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The response of two-dimensional electron gas to temperature gradient in perpendicular magnetic 
field under steady-state microwave irradiation is studied theoretically. The electric currents in¬ 
duced by temperature gradient and the thermopower coefficients are calculated taking into account 
both diffusive and phonon-drag mechanisms. The modification of thermopower by microwaves takes 
place because of Landau quantization of electron energy spectrum and is governed by the microscopic 
mechanisms which are similar to those responsible for microwave-induced oscillations of electrical 
resistivity. The magnetic-field dependence of microwave-induced corrections to phonon-drag ther¬ 
mopower is determined by mixing of phonon resonance frequencies with radiation frequency, which 
leads to interference oscillations. The transverse thermopower is modified by microwave irradiation 
much stronger than the longitudinal one. Apart from showing prominent microwave-induced oscil¬ 
lations as a function of magnetic field, the transverse thermopower appears to be highly sensitive 
to the direction of linear polarization of microwave radiation. 
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I. INTRODUCTION 

Electron transport in two-dimensional (2D) electron 
systems placed in a perpendicular magnetic field remains 
one of the most important subjects in the condensed 
matter physics. Recently, it was established that a va¬ 
riety of interesting transport phenomena takes place^ in 
the region of moderately strong magnetic field, where 
the Shubnikov-de Haas oscillations of the electrical re¬ 
sistivity are suppressed because of thermal smearing of 
the Fermi level. In particular, there are several kinds of 
magnetoresistance oscillations^ observed in high-mobility 
2D systems such as GaAs quantum wells, strained Ge 
quantum wells, and electrons on liquid helium surface. 
Among these phenomena, the microwave-induced resis¬ 
tance oscillations (MIRO) appearing under microwave 
(MW) irradiation of 2D electron gas^“^ are studied most 
extensively. Their origin is briefly described as follows. 
In the presence of the MW excitation, when absorption 
and emission of radiation quanta by electron system take 
place, both the distribution function and scattering prob¬ 
abilities of electrons are modified. These modifications 
correlate with the oscillating density of electron states 
owing to Landau quantization in the magnetic field B, 
thus leading to corresponding oscillating contributions to 
resistivity determined by the ratio of the radiation fre¬ 
quency w to the cyclotron frequency ujc = \e\B/mc. The 
period and phase of MIRO, as well as the temperature 
and power dependence of their amplitudes are in agree¬ 
ment with this physical picture supported by a detailed 
consideration of microscopic mechanisms of MIRO in the 
past years.According to both experiment and the¬ 
ory, the MW irradiation strongly affects the longitudinal 
(dissipative) resistivity and has a much weaker effect on 
the transverse (Hall) resistivity. More recent experiments 
uncover the existence of small corrections, sensitive to 
the direction of electric field of microwaves (MW polar¬ 


ization), to both longitudinal and Hall resistivities, 
also in general agreement with the theory. 

Apart from its influence on electrical resistivity, the 
MW excitation is expected to modify other transport co¬ 
efficients of 2D electrons, for the same reasons as ex¬ 
plained above. The thermoelectric coefficients are of 
special interest in this connection. The study of ther¬ 
moelectric phenomena in magnetic fields has a long his¬ 
tory, and the fundamentals of this topic, with appli¬ 
cations to bulk conductors, are reviewed in Ref. 14. 
The electrical response to temperature gradient VT 
is described by the longitudinal (Seebeck) and trans¬ 
verse (Nernst-Ettingshausen) components of thermoelec¬ 
tric power (briefly, thermopower). These coefficients are 
determined by two mechanisms: the diffusive one, when 
electrons are directly driven by the diffusion force due to 
temperature gradient in electron gas, and phonon drag 
one, when electrons are driven by a frictional force be¬ 
tween them and phonons propagating along the tempera¬ 
ture gradient. Contribution of both these mechanisms in 
magnetothermopower of 2D electron systems has been 
studied in a number of theoretical and experimental 
works^®“^^ (see also review paper Ref. 22 and references 
therein). The quantum effects are commonly observed in 
strong magnetic fields, where Shubnikov-de Haas oscilla¬ 
tions of thermopower coefficients take place.In high- 
mobility GaAs quantum wells, the phonon-drag ther¬ 
mopower shows another kind of quantum oscillations, re¬ 
lated to resonant phonon-assisted scattering of electrons 
between Landau levels.This occurs in the region of 
moderately strong magnetic fields, below the onset of the 
Shubnikov-de Haas oscillations, which is favorable for ob¬ 
servation of MW-induced quantum effects. 

There are two main ways in which the MW irradiation 
can influence the thermopower. First, this irradiation 
leads to non-equilibrium electron distribution that has 
non-trivial dependence not only on electron energy but 


2 


also on temperature of electron gas. Both the diffusive 
and phonon-drag contributions to thermoelectric coeffi¬ 
cients should be sensitive to such changes. Next, the 
MW irradiation in the presence of magnetic field consid¬ 
erably influences electron-phonon scattering. This causes 
an effect on electrical resistivity^^ under conditions when 
probability of electron-phonon scattering is comparable 
to that of elastic scattering by impurities. At the tem¬ 
peratures below 4.2 K these conditions are realized only 
in very pure 2D systems. In contrast, the effect of mi¬ 
crowaves on electron-phonon scattering is always impor¬ 
tant for thermoelectric properties, since the phonon drag 
mechanism determined by this scattering gives a very 
significant,^^ if not a major, contribution to thermopower 
of 2D electrons. 

The above consideration also suggests that in spite of 
the same microscopic mechanisms involved in both cases, 
the effect of microwaves on magnetothermoelectric coef¬ 
ficients of 2D electrons should be different from their ef¬ 
fect on magnetoresistance. The classical Mott relation 
between the diffusive current responses to temperature 
gradient and to electric field is not expected to be valid 
under MW excitation, even for the moderately strong 
magnetic fields. Moreover, one may presume that both 
longitudinal and transverse components of thermopower 
oscillate with magnetic field in the way different from 
MIRO, and their dependence on MW polarization is also 
different. Therefore, there is enough motivation for a the¬ 
oretical study of the influence of MW irradiation on ther¬ 
mopower of 2D electron systems in perpendicular mag¬ 
netic field. The present paper is devoted to this previ¬ 
ously unexplored problem. 

In the linear response regime considered in the follow¬ 
ing, the current density j is given by the general expres¬ 
sion 

j = dE - ^vr, (1) 

where E is the electric field in the plane (x, y) of the 
2D electron system. It is assumed that the 2D system 
is macroscopically homogeneous so that the chemical po¬ 
tential y does not depend on 2D coordinate r. Under 
condition when no conduction currents flow in the sys¬ 
tem, one gets E = dVT. The thermopower tensor a 
describes the voltage drop as a result of temperature gra¬ 
dient. It is given by d = p/3, where the resistivity ten¬ 
sor p is the matrix inverse of the conductivity tensor a. 
The theoretical approach presented below is based upon 
calculation of thermoelectric tensor /3 in the presence of 
ac field of microwaves by using the method of quantum 
Boltzmann equation^’®do.23 established in the previous 
calculations of the conductivity tensor a. As both /3 
and a are known, the thermopower is found straightfor¬ 
wardly. The results are presented for the case of moder¬ 
ately strong magnetic field, when the Shubnikov-de Haas 
oscillations are still suppressed, but quantum oscillations 
due to Landau quantization exist in high-mobility 2D 
systems. Such oscillations are caused by inelastic scat¬ 
tering of electrons between Landau levels as a result of 


electron interaction with acoustic phonons of a resonant 
frequency ojph (magnetophonon effect^^’^"*”^®) and with 
microwaves of frequency w. These two kinds of inelastic 
processes actually interfere, leading to combined frequen¬ 
cies ujph ± w whose ratio to Wc determines the periodicity 
of the quantum magnetooscillations.^^ As shown below, 
such oscillations exist in both longitudinal and transverse 
thermopower caused by the phonon drag, while the diffu¬ 
sion part of the thermopower follows the MIRO period¬ 
icity determined by the single frequency uj. The phonon- 
drag part of MW-induced contribution to transverse ther¬ 
mopower is found to be comparable with that of longi¬ 
tudinal thermopower. Since the transverse thermopower 
is much smaller than the longitudinal one in classically 
strong magnetic fields, it is dramatically affected by MW 
irradiation, demonstrating giant microwave-induced os¬ 
cillations and a high sensitivity to the direction of MW 
polarization. 

The paper is organized as follows. Section II de¬ 
scribes the main formalism including description of ac 
electric field generated by incident electromagnetic ra¬ 
diation, electric current in the presence of temperature 
gradient, and kinetic equation for 2D electrons with colli¬ 
sion integrals for electron-impurity and electron-phonon 
interactions. In Section III the kinetic equation is solved 
and the tensor jS is presented and discussed both for the 
equilibrium case and under MW irradiation. Section IV 
contains expressions for longitudinal and transverse com¬ 
ponents of thermopower tensor d, their discussion, and 
presentation of the results of numerical calculations of 
these components as functions of magnetic field and po¬ 
larization angle. More discussion and concluding remarks 
are given in the last section. 


II. GENERAL FORMALISM 

Throughout the paper, one uses the system of units 
where Planck’s constant h and Boltzmann constant ks 
are both set to unity. The electron spectrum is assumed 
to be isotropic and parabolic, with effective mass m. The 
Zeeman splitting of electron states is neglected. 

Consider a monochromatic electromagnetic wave nor¬ 
mally incident on the surface containing a 2D layer (the 
direction of incidence coincides with the direction of the 
magnetic field, along the z axis). The electric field of this 
wave near the layer is written, in the general form, as 





where e is the polarization vector. The second part of 
this equation represents the wave as a sum of two circu¬ 
larly polarized waves, e± = (e^, ± iey)l'/2 = x 

is the angle between the main axis of polarization of E^*^ 
and the x axis, and k± are real numbers characterizing 
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ellipticity of the incident wave (they are normalized ac¬ 
cording to = 1). A circular polarization means 

that either or is equal to zero. In the case of linear 
polarization, k+ = K- = so that e± = 

The screening of electromagnetic wave due to the pres¬ 
ence of free carriers in the 2D layer changes the polariza¬ 
tion angle and ellipticity,so the electric field in the 

U\ 

layer, Ej, differs from EJ and has the following form: 


Et = :^Re 




-|-(w -|- Wc)s-|- 


1 

—i 



( 3 ) 


where 


s± 


e± 

LO ±UIc + itOp 


( 4 ) 


Here cOp is the radiative decay rate that determines 
the cyclotron line broadening because of electrodynamic 
screening effect. It is given by uip = 2Tre'^nsjmc'/^, 
where is the electron density, -x/e* = (1-1- 
and e is the dielectric permittivity of the sample ma¬ 
terial. Next, In Eqs. (3) and (4), it 

is assumed that transport relaxation rate, which de¬ 
termines electron mobility, is much smaller than either 
|w ± Wcl or iOp. The relation vtr <C Wp is a very good ap¬ 
proximation for high-mobility samples with typical elec¬ 
tron density Us > 10^^ cm“^. 

Apart from the ac field Et, the electron system is 
driven by a weak static (dc) held E. To take into ac¬ 
count the inhuence of both these helds on 2D electrons, 
it is very convenient to use a transition to the moving 
coordinate frame (see Ref. 10 and references therein). 
Then, the quantum kinetic equation for electrons derived 
by using Keldysh formalism for non-equilibrium electron 
systems (see details in Refs. 10 and 23) contains the ef¬ 
fect of external helds only in the collision integral. The 
radiation power is assumed to be weak enough to neglect 
the inhuence of microwaves on the energy spectrum of 
electrons: the spectrum remains isotropic and the den¬ 
sity of states is not affected by the radiation. Further, the 
magnetic held is assumed to be weak enough so there is 
a large number of Landau levels under the Fermi energy. 
The kinetic equation written for the electron distribu¬ 
tion function fg^p averaged over the period 27r/a; takes 
the form 




V feip + 


dfeif 
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where e is the electron energy, = ^^(cos v?, sin (/?) 
with pg = \/2me is the electron momentum in the 2D 
layer plane, and ip is the angle of this momentum. Since 
the dependence of all quantities on the spatial coordi¬ 
nate r is considered as a parametric one, the coordinate 
index at the distribution function and collision integrals 
is omitted. The density of electric current is given by the 


expression 

i=^JdeDgJ ^ Pecp/ ecp - cr±eE - ceVMz, (6) 

where cr± = e^UsImuJc = \e\nsc/B is the classical Hall 
conductivity and Dg is the density of electron states ex¬ 
pressed in the units to/tt. Next, e = ^ ^ ^ is the 

antisymmetric unit matrix in the space of 2D Cartesian 
indices. The last term in the expression (6) is given by 
the spatial gradient of magnetic moment M of electrons 
per unit square. This moment arises because of diamag¬ 
netic currents circulating in the electron system. Actu¬ 
ally, the last term in Eq. (6) does not contribute to the 
total current across any hnite sample. However, the ne¬ 
cessity of taking into account this term (its bulk analogue 
is —c[V X M]) in the expression for the local current den¬ 
sity has been emphasized in studies of magnetothermo¬ 
electric phenomena a long time ago.^^’^^ Being expressed 
through the distribution function, the magnetic moment 
comprises two terms: 

Mg = -^ j ds[Dge-ng]fg, (7) 

where fg is the isotropic (averaged over ip) part of 
electron distribution function, and Hg = ds'Dgi 
is the antiderivative of Dg. In the ideal 2D electron 
system, the first and the second terms in Mz corre¬ 
spond to magnetization due to bulk and edge currents, 
respectively.^^ In the case of the equilibrium Fermi dis¬ 
tribution function fg = {exp[(£ — p)/T] + 1}“^, it 
is easy to transform Eq. (7) to a well-known ther¬ 
modynamic expression Mz = —dfl/dB, where D = 
— {Tmlir) f deDg\o.{l-\-eyip[{p — e)/T]} is the thermo¬ 
dynamic potential per unit area. 

In the absence of any collisions, Jgp = 0, the local 
current is non-dissipative, j = j*'"\ where 

=-^f dsIlgeVfg - UicE. (8) 

If coordinate dependence of fg exists solely due to tem¬ 
perature gradient, one has Vfg = {dfg/dT)VT. The 
integral term in Eq. (8) is reduced to non-dissipative 
thermoelectric current —fi'S/T flowing perpendicular to 
VT, with P = {cmlTrB) f deTIg(dfg/dT)i. If chemical 
potential p entering fg also depends on coordinate, the 
integral in Eq. (8) produces an additional term propor¬ 
tional to V/i. This term, with the aid of the identity 
dfg/dp = —dfg/de, can be combined with the last term 
of Eq. (8), leading to the form a±eV(, where f = ^+p/e 
is the electrochemical potential and $ is the electrostatic 
potential determining the electric field E = — Vd". The 
electric field or, in general, the gradient of electrochem¬ 
ical potential induced as a result of a temperature gra¬ 
dient is derived from the expression = 0. This leads 
to diagonal thermopower tensor d = fa, where 1 is the 
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unit 2 x 2 matrix and 


m 

a = - 

ireus 


deU, 


dT' 


(9) 


Substituting the equilibrium distribution function into 
Eq. (9), one gets a well-known result 


S 

\e\ns’ 



( 10 ) 


where S is the entropy of 2D electron gas per unit area. 
For strongly degenerate electron gas, jjL = ep, one has 
S = {'JT^/3)nsT/eF- 

The collision integrals and standing in Eq. 
(5) describe, respectively, electron-impurity and electron- 
phonon scattering:^^ 


/•27r J,./ °° 

^ E K|qen|)[Jn(|R.-qe„|)]^ 

n—— oo 

X -^£+na;+7n [fe-\-nu}-\-'yn^' fscp], (11) 


jph _ 
'^Sip ~ 
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dip' ^ /■“ dq, ^ 

27r E / 9^"^ E 


Jo ^ t/—oo _ 

• qe„|)]^[(.^AQ- + fsip) 

^ f €—uj^Q — +nuj+^P p' (-^AQ “b Ofsp] 

^ ^s-LO^^-+nLO+'rP + -^AQ+["^ndR-w ’ q)!dl)] 
^[(■^A-Q+ + 1 - fscp)fspi^^^_^_pnuj+'f+ip' 

— Nx-Q+feip]Dgp^^^^pnuj+'f+'\^ ( 12 ) 


where J„ is the Bessel function, ^{q) = mw(q) is the 
isotropic elastic scattering rate expressed through the 
Fourier transform w{q) of the correlation function of ran¬ 
dom potential of impurities, q^n = Pep — Ps+nujp' is the 
momentum transferred in scattering in the presence of ac 
field, and R^j is a complex vector describing coupling of 
electron system to this field: 


R. = -^(s+ + s-,(s+-s_)A). (13) 

V Zmoj 

The interaction with phonons is considered under ap¬ 
proximation of bulk phonon modes. The phonons are 
characterized by the mode index A and three-dimensional 
phonon momentum Q with out-of-plane component qz- 
The squared matrix element of electron-phonon interac¬ 
tion potential is represented as Mxq = CxQlq^- The 
squared overlap integral Iq^ = \ ( 0 |e*®^^| 0 ) p is deter¬ 
mined by the confinement potential which defines the 
ground state of 2D electrons, |0). The function Cxq char¬ 
acterizes electron-phonon scattering in the bulk. The 
in-plane momenta transferred in electron-phonon colli¬ 
sions, q^„, are found from the equation q^„ = — 

Pe±u,^^±+nuip', where = (q^^,^^) and waq is the 
phonon frequency. The effect of the static electric field 


on the collision integrals is given by the energies 7 „ = 
Yo-qen and 7 ^ = VD-q^„, where Yp = c[Ex B]/^^ = 
{c/B)eEi is the drift velocity in the crossed electric and 
magnetic fields. 

In the case of electrons interacting with long- 
wavelength acoustic phonons in cubic lattice, the expres¬ 
sions for Cxq and dynamical equations needed for deter¬ 
mination of ijJxq are the following: 


Caq = 


1 


2 PmWaq 


p2 E eAQiSAQi?! 


Wi 


^ (ehi4)^ ^ 

“r ^ j'k'^xqk^xqk'qiqjqpqj' 

V 


ijk,i' j'k' 


, ( 14 ) 


E [^b(Q) - exQj = 0 , (15) 

3 

^ iji .03 — [(^11 ^44)0’^ “b 044!^ 

+ ici2 + C44)qiqj{l-6ij). (16) 


Here V is the deformation potential constant, /114 is the 
piezoelectric coupling constant, and Pm is the material 
density. The sums are taken over Cartesian coordinate 
indices. The coefficient Kijk is equal to unity if all the 
indices fc are different and equal to zero otherwise. 
Next, e^Qi are the components of the unit vector of the 
mode polarization, and KijlQ) is the dynamical matrix 
expressed through the elastic constants cn, C 12 and C 44 . 

Finally, Nxq in Eq. (12) is the distribution function of 
phonons. In the presence of thermal gradients this func¬ 
tion depends not only on the frequency ujxq, but also 
on the direction of Q. In the general case, Nxq can be 
represented as a sum of symmetric (s) and antisymmet¬ 
ric (a) parts satisfying the relations = N^q and 

iV^_Q = —N^q, respectively. The drag of electrons by 
phonons is caused by the antisymmetric part. In the lin¬ 
ear regime, iV“ is proportional to VT while is reduced 
to the equilibrium distribution function. In particular, 
one often uses the following form:^^ 


Nxq = N^^^ 


^Nuxq ^XQ 
dujxQ T 


taUaq • VT, 


(17) 


obtained from a linearized kinetic equation for phonons 
in the relaxation time approximation. Here N^^^^ = 
[exp(a;AQ/T) — 1]“^ is the equilibrium (Planck) distri¬ 
bution function, ta is the relaxation time of phonons, 
and uaq = Swaq/SQ is the phonon group velocity. No¬ 
tice that a simple expression Uaq = sxQ/Q relating the 
group velocity to the sound velocity sx is valid only in 
the isotropic approximation. For elastic waves in real cu¬ 
bic crystals the direction of uaq is not generally coincide 
with the direction of Q, though the symmetry relation 
ua-q = —Uaq is always valid. Substituting Eq. (17) 
into the expression for the collision integral it is 
convenient to write the latter as a sum of two parts, 



( 18 ) 
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where contains the equlibrium phonon distribution 

iV^^p only, while is determined by the anisotropic 

non-equilibrium correction to phonon distribution [sec¬ 
ond term in Eq. (17)] and is proportional to VT. The 
second term in Eq. (18) is responsible for the phonon 
drag contribution to electric current. 

By using Eqs. (11) and (12), one can directly check 
the identity f deD^ f = 0 expressing electron con¬ 

servation requirement. It is worth to emphasize that the 
collision integrals Eqs. (11) and (12) are written in the 
general form valid for arbitrary relation between radi¬ 
ation frequency w, phonon frequency waq and electron 
energy e. In Ref. 23 the collision integrals are written 
in a simpler form valid under the assumptions a; <C e 
and waq «C £. For degenerate electron gas, the electrons 
contributing to electric current have energies £ close to 
the Fermi energy, and these assumptions usually work 
very good for microwave frequencies and acoustic phonon 
scattering. However, in the problem of diffusive ther¬ 
mocurrent an extra accuracy is required, so the terms of 
the first order in w/£ are to be retained at least in the 
isotropic part of electron-impurity collision integral (see 
Eq. (25) below). 


III. SOLUTION OF KINETIC EQUATION 


When searching for the response to temperature gradi¬ 
ents only, the effect of the dc field in the collision integrals 
is omitted, 7 n = 7 )^ = 0. It is also assumed that the main 
cause of momentum relaxation of electrons comes from 
electron-impurity scattering rather than from electron- 
phonon scattering. In GaAs quantum wells with electron 
mobility of about 10® cm^/V s this approximation holds 
al low temperatures T < 10 K (for GaAs quantum well 
of typical width 20 nm the phonon-limited mobility is es¬ 
timated as 1.3 X 10^ cm^/V s at T = 4.2 K and 3.8 x 10® 
cm^/V s at T = 10 K). Thus, one may neglect the contri¬ 
bution in comparison to T™, but the contribution 

leading to phonon drag must be retained. It is 
convenient to expand the distribution function in the an¬ 
gular harmonics, The electric current 

density given by Eq. (6) is determined by the compo¬ 
nents with k = ±1. Only the effects linear in MW power 
are considered below. The distribution function is repre¬ 
sented as a sum of two terms, where 

is proportional to MW power. For fc 0 the first term 
is given by the expression comprising the diffusive and 
phonon-drag parts: 


AO) _ _ 

sk 
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ikujc + \ 2m 


Pe (9fek+l 


dT 


V+T 


dfek- 

dT 


-V_T -b 


c 2 'k 


^ gi(fe' 


-k)ip 


C.M I ^ - /efc')||, (19) 


where V± = ± zVy, = V0[1 — cos(fc6l)] (the line 

over the expression denotes angular averaging), and vg = 
jz[2p£ sin(0/2)]. The integral operator M is proportional 
to temperature gradient and defined as 






( 20 ) 


with F{x) = [a;/sinh(a:)]^. It is taken into account that 
waq <C £, which allows one to use the quasielastic approx¬ 
imation, when the transferred 2D momentum is re¬ 
placed by q, with absolute value q = 2ps sin(6l/2) depend¬ 
ing on electron energy and scattering angle 9 = (p — ip'. 
The angle of the vector q is = 7r/2 -b <(), where 
(j) = {p + p')/2. The phonon frequency can be written as 
waq = sxqQ, where Q = \/+ ql and saq is the sound 
velocity that depends on the mode index and direction 
of vector Q. If the quantum well is grown in the [001] 
crystallographic direction, as assumed in the following, 
both waq and M\q are periodic in pq with the period 
7r/2. 

To find with the accuracy up to the linear terms 

in MW power, only the contributions with low-order, 
|n| < 1, Bessel functions Jn are to be taken in Eqs. (11) 
and (12). Physically, this corresponds to a neglect of 
multi-photon absorption processes. If fc 0, then 


r{MW)_ 7A(£)/4 

~ ikujc + i'^’^W, 


f^-dp f^^dp' 


x[l-6e2*'^- 

k' n=±l I 


ve 


x[D,+^Ms+nu,k'e-''^’^- f,k') - D,f,k'{e-'^'^- 1 )] 
~-M •! ^ ] I ^Di:-li^^^+nuiifs-luJxQ+nujk'e- ~ fek') 

- U')]}} , ( 21 ) 


—ik' 9 






where 

P.(£) = ^^(|s+p + ls-n (22) 

is the dimensionless function proportional to MW power 
(see Eq. (4) for definition of s±) and 


b = 


S_Sj 


(|s+l 


S_P) 


(23) 


is a complex dimensionless coefficient which depends on 
the direction of MW polarization and determines the sen¬ 
sitivity of transport properties of electrons to this di¬ 
rection. The neglect of multi-photon processes implies 
Pujis) "C 1. The expression (21) comprises both electron- 
impurity and electron-phonon parts, though only the 
electron-phonon part is essential below. 
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To find the isotropic (fc = 0) part of the distribution 
function, it is necessary to include electron-electron scat¬ 
tering into consideration. Though the corresponding col¬ 
lision integral is not written in Eq. (5) explicitly, it 
can be found in Refs. 9 and 23. The kinetic equation is 
written as 


J^rn ^ jph(0) jee ^ 


(24) 


where only the isotropic part of the distribution func¬ 
tion is retained under the collision integrals. It is essen¬ 
tial that J®® is not affected by MW irradiation, while 
J*™ is non-zero only in the presence of MW irradia¬ 
tion. The distribution is represented as a sum of 
smooth part and rapidly oscillating part . The 

smooth part is controlled by electron-electron scattering 
and, therefore, can be approximated by a heated Fermi 
distribution with effective electron temperature Tg, the 
latter is to be found from the energy balance equation 
/ deDi.e[Jl"^ + = 0. For oscillating part, one gets 

the following expression: 


y; (l + ^(1 - Z„)) 

n=±l 

x6D,+„Me+nu - fe), (25) 

where = 0^—1, vtr = is the transport 

relaxation rate, and 


2tr = 


d In Ttr 

91n£ 


(26) 


is the logarithmic derivative of the transport time over 
energy. The inelastic scattering time Tin entering Eq. 
(25) describes relaxation of the isotropic oscillating part 
of electron distribution.® This relaxation is caused mostly 
by electron-electron scattering and scales with tempera¬ 
ture as Tin oc T~'^. 

The electric current is given by Eq. (6), where the 
distribution function, found from Eqs. (19), (21), and 
(25) with the accuracy up to the terms linear in both 
Vr and P^, is substituted. The thermoelectric tensor /3 
determining the thermocurrent is represented below as a 
sum of four parts: 


/3 = , (27) 

where diffusive (d) and phonon drag (p) parts are writ¬ 
ten separately. Two first terms correspond to dark ther¬ 
mocurrent, in the absence of MW irradiation, while the 
next two terms are MW-induced corrections. While 
and are determined only by from Eq. (19), the 
MW-induced parts are found in a more elaborate way, by 
combining together the results given by Eqs. (19), (21), 
and (25), as described in Subsection B. 


A. Dark thermocurrent 

In the absence of MWs, the linear response to temper¬ 
ature gradient is found from Eq. (19) for k = ±1 with 


isotropic (fc' = 0) distribution functions substituted in 
the right-hand side. The thermoelectric coefficients are 
given by the following expressions: 


2(0) |e| f dfP 

~ T J dT, ul + vlDl 


2(0) _ M /w UJcSDe^- 

2^7 ^ ujl + ulDl 

xA I ^ - /i°7o^AQ| , 


(28) 


(29) 


where Vn is the integral operator defined as 



X (1 - cos OYwPMxqtxF (^) (30) 

The matrices given by Eqs. (28) and (29) contain diag¬ 
onal symmetric (oc 1) and non-diagonal antisymmetric 
(oc e) parts, so their symmetry is the same as the sym¬ 
metry of the electrical conductivity. 

The expressions (28) and (29) describe the thermoelec¬ 
tric tensor in a wide region of temperatures and mag¬ 
netic fields. Quantum osicllations of P^^ and occur 
because of the oscillating dependence of the density of 
states, Dg. In the following, the approximation of over¬ 
lapping Landau levels is used: = 1 — 2dcos{2Tre/uJc), 

where d = exp(— tt/IwcIt) is the Dingle factor {d <C 1) 
and T is the quantum lifetime of electrons, given at 
low temperatures by r = l/l7e- Apart from the condi¬ 
tion d <C 1, the validity of the expression for im¬ 
plies ST ^ 1. Under the same requirements, fig = 
e—(uJc/Tr)dsm(2Tre/uJc)- The integrals over energy in Eqs. 
(28) and (29) are calculated below under the assumption 
of strongly degenerate electron gas, and the quantum ef¬ 
fects up to the second order in the Dingle factors are 
retained. To take into account energy dependence of the 
Dingle factor due to energy dependence of r, the loga¬ 
rithmic derivative Z = 91nr/i91ne is introduced. The 
diffusive part is given by the following expression: 


= 


5|Te 


-l-6d cos 


3(a;2 

2TreF B 


— 12d-^^Bsin 
nTe 


OJc X 

2TTeF 




e - Vtr 


tOc 


1 + 


2Ztrvh 




+ ^tr 


l-Zi 


OJc 


2dM 1 - Ztr + 


2 2 
Wc - vfn 
2 , 2 
+ I'tr 

2ttZ 


\0Jcp 


1 (31) 


with X = 2'tPTcIoJc and B = 9(A/sinhA)/9A = 
(1 — A cothA)/sinhA. All energy-dependent quanti¬ 
ties, namely Vic-, t, Zt^, and Z, in Eq. (31) are taken 
&t e = Ef- The classical terms and the quantum term 
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proportional to d in the diagonal part of have been 
reported previously.^^ 

Calculating the phonon-drag part from Eq. (29) under 
the same approximations, one gets the result 



-i^tr [Ti{l + 2d^)+4d'^T^i] i 


-ddEsi . ^ y cos ^ojce - (3/2)r'tri] 1 • (32) 

Sinn Ji. cjc J 


Similar to Eq. (31), this expression contains both classi¬ 
cal terms and quantum terms proportional to d and d^. 
The dimensionless functions E^ used here and below are 
defined as 



= V 

' n 


cos 


1 

2TruixQ 


2-n-uJxQ 


sm 


2TruJxQ 


(33) 


where denotes Vn e = ep- The function Ti deter¬ 
mines the classical contribution^^ to phonon-drag ther¬ 
moelectric response and does not depend on the mag¬ 
netic field. This contribution has been considered pre¬ 
viously in the isotropic approximation for phonon spec¬ 
trum, when there are one longitudinal phonon branch 
with velocity si and two transverse branches with ve¬ 
locity St- For high temperatures, when T exceeds both 
s\Pf and nsx/a {pp is the Fermi momentum and a is 
the quantum well width), Ti is temperature-independent. 
At low temperatures, T <C s\pp, the Bloch-Gruneisen 
transport regime is realized, when electron scattering by 
phonons occurs at small angles, 0^1. In this regime, 
Ti scales with temperature as (or as if only the 
deformation-potential mechanism of electron-phonon in¬ 
teraction is present). The functions Ten and T^n standing 
in the quantum contributions depend on the magnetic 
field and can be analytically calculated only in certain 
limits (see Appendix). 

The terms oc d in Eqs. (31) and (32) describe the 
Shubnikov-de Haas oscillations of the thermocurrent. 
The forthcoming consideration, however, is focused at 
the case of |wc| 27r^Te, which means that A/sinhX 
is exponentially small so the Shubnikov-de Haas oscil¬ 
lations are suppressed and the quantum corrections are 
given by the terms oc d^ only. Using Ug = mep/p, one 
may check that the tensor (31) under these conditions 
satisfies the Mott relation ,5^°^ = —(p'^Te/i\e\){da/dep), 
where d = iod — cup is the conductivity tensor whose 
components are Ud = e^nsVtr{l + 2d^)/[m{uj1 -I- and 
Up = + I'tr)]- The quantum corrections 

both in these expressions and in Eq. (31) are essential 
only in the classically strong magnetic fields, so the terms 
oc {vtr/‘^c)d^ are neglected in comparison to the terms 
oc d^. 

The diffusive thermoelectric coefficients do not oscil¬ 
late before the onset of Shubnikov-de Haas oscillations. 


In contrast, quantum magnetooscillations of phonon- 
drag thermoelectric coefficients persist under the as¬ 
sumed condition |wc| «C 27r^Te, because of the presence of 
Tci in . Indeed, the oscillating nature of the function 
cos{2pujxci/uJc) is not completely washed out after the 
integration under V. The major contribution to such in¬ 
tegrals comes from the region of variables around qz = 0 
and 0 = TT, which physically corresponds to backscat- 
tering of electrons as a result of emission or absorp¬ 
tion of phonons moving in the quantum well plane, the 
wavenumber of these phonons is close to 2pp. Thus, there 
exist resonant phonon frequencies, roughly estimated as 
“^PfSx, which lead to magnetooscillations of phonon- 
drag thermopower observed^*^ in high-mobility samples. 
With decreasing temperature, the oscillations are ex¬ 
ponentially suppressed in the Bloch-Gruneisen regime 
(see Appendix). The same kind of oscillations is ob¬ 
served in electrical resistivity, they are known as acous¬ 
tic magnetophonon oscillations or phonon-induced resis¬ 
tance oscillations. 


B. Microwave-induced thermocurrent 

The distribution functions f^i and fe-\ determining 
the electric current under MW irradiation are to be found 
up to the terms linear in Pai{s). There are two types of 
such MW-induced contributions. The direct ones are ob¬ 
tained in two ways: i) by calculating from Eq. 

(21), where the isotropic distribution function is re¬ 
tained under the integral (only the phonon part is essen¬ 
tial), and ii) by calculating from Eq. (19), where 

the isotropic MW-induced distribution function 
is placed in the right-hand side. The indirect contribu¬ 
tions assume calculation of and by substi¬ 
tuting anisotropic parts of and , respectively, 

in the right-hand sides of Eq. (19) and Eq. (21). Sim¬ 
ilar technique has been used for calculation of the MW- 
induced conductivity. Following the notations of Ref. 
10 , one may denote the direct contributions (i) and (ii) 
as the ’’displacement” and ’’inelastic” ones, respectively, 
and the indirect contributions as the ’’quadrupole” ones. 
Strictly speaking, there exists one more indirect contribu¬ 
tion called the ’’photovoltaic” one,^° which is determined 
by the MW-generated time-dependent part of the distri¬ 
bution function and cannot be obtained from the kinetic 
equation Eq. (5) because the latter is written for time- 
independent fe^p. The indirect contributions to j5 begin 
with the terms of the order vtrjcjc compared to direct 
contributions. Since all the MW-induced contributions 
are of quantum nature and important only in the region 
of classically strong magnetic field, ujc ^ vtr, the indirect 
contributions are less significant than the direct ones and 
can be safely neglected in the thermopower coefficients 
presented in the next section. Therefore, the attention 
below is focused at the direct contributions only. 



The current is calculated in the regime when electron 
gas is degenerate. Within the required accuracy, the so¬ 
lution of Eq. (25) is given by the following expression: 


jiMW) ^ 

Z CJc 


X (/i+l - /i-l) - cos 


27r£ . 27ra; 

sin- 

OJc 

2ttui 




X 

27rZ 

|Wc|t 


2716 

p(0) 


• COS ■ 


UJc 


f(0) _|_ f(0) _ 9f(0) 

Je-\-uj ' Je—uj 


(AO) _ AO) ( 
\J s-\-uj Je—Lj) 


(34) 


The first term of this expression gives the main con¬ 
tribution sufficient for calculation of the MW-induced 
resistance.® The second term represents a correction of 
the order w/e, which is necessary for calculation of the 
diffusive thermopower. The term proportional to the 
factor sin(27rw/a;c) in Eq. (34) also enters Eqs. (35), 
(37), (43), (45), and (48) below, where the contribu¬ 
tion of inelastic mechanism is present. This factor re¬ 
flects the property® that the strongest modification of 
the electron distribution function under MW irradiation 
in the presence of weak Landau quantization occurs when 
oj/ lOc = n ± 1/4 (n is an integer). The correction pro¬ 
portional to the factor cos(27rw/a;c) appears because the 
resonance absorption of MW radiation at ujjojc = n also 
has an effect on the distribution function. 

The consideration below assumes the approximation 
\ujc\ «C 27r®Te, when Shubnikov-de Haas oscillations are 
thermally averaged out. This dramatically simplifies cal¬ 
culation of the integrals over energy because one can av¬ 
erage the products of rapidly oscillating functions such 
as and over the period Wc before integration 

over the energy. After substituting Eq. (34) into the first 
part of the right-hand side of Eq. (19) and calculating 
the current according to Eq. (6) (one may equally use 
Eq. (28) with replaced by /i^'^^), the diffusive part 
of /3 takes the form 


X 


AMW) 

Pd 


\e\(fnnVtrOj‘^Pa>%n 
nTeiuj'^ + 


. 27TUJ ^ 

-Sin-e — 

271LO OJc 


^ \ 27T0J- 

^tr) COS 1 
OJc 


(35) 


where all energy-dependent quantities are taken at e = 
£f, in particular, = PujOf)- The main contribution 
to the derivative over temperature in Eq. (19) comes 
from temperature dependence of the inelastic relaxation 
time, expressed through the logarithmic derivative 




51nTi„ 

ainTe “ 


(36) 


The factor sin(27ru;/wc) typical for MW-induced 
conductivity® does not appear in the diagonal part of 
thermoelectric tensor Eq. (35), because of different de¬ 
pendence of the diffusive thermoelectric current on elec¬ 
tron energy distribution as compared to the drift current. 


The first term of is averaged out in the diagonal 

components of , while the second term of ^ 

proportional to cos{2Truj/uJc), survives this averaging. 

For phonon-drag part of /3 the result is the following: 


pi^MW) 


where 


2|e|ngd^P^ 
+ Vtr) 


27ru; , 27rw 

VtrTint si- Sin- 

Uc OJc 

X [—OJci -\- {i/2')l'tr^\ 

. -rs ■ Z J^OJ ^ FOJ . 2'K0j\ 

-|-(rc 2 Sin-|-rs 2 —sm- 

OJc OJc OJc ) 

X [wc(—e -l- go) + 2t'tr(i + ^o)] 

1 X, .9 T^OJ ~ TTOJ . 2Troj\ 

-I- rc 2 sm-h rs 2 -- 1 


OJc 


sm ■ 

OJc OJc ) 


X [ojcgi + 2Ftrhi]j , (37) 


go = b'ax + b"a^, gi = b'a^ - P'&z, 

ho = b'az - b"a3,, hi = b'az A b''(Jx, (38) 




1 0 
0 -1 


and Ux 


0 1 
1 0 


are the Pauli matri¬ 


ces, while b' and b" denote real and imaginary parts of b 
[see Eq. (23)], respectively. The quantities Ti differ from 
Ti by placing the factor 


cos(4(pg) 


sin(4(p5) 



4p| sim 0/2 


— ^(39) 

OJ\Q Oipq 


under the integral operator P in Eq. (33). In the general 
case, both terms in Eq. (39) are essential for calcula¬ 
tion of Tj. In the isotropic approximation for phonon 
spectrum the second term in Eq. (39) vanishes, but 
is still nonzero, because the piezoelectric-potential part 
of M\q remains angular-dependent. If the anisotropy of 
phonon spectrum is weak, the second term in Eq. (39) 
can be neglected in calculation of piezoelectric potential 
contribution. 

The expression (37) includes contributions from both 
inelastic (first term) and displacement (second and third 
terms) mechanisms. The inelastic-mechanism contribu¬ 
tion can be obtained from Eq. (29) after replacing 
by j/^^) ^ displacement-mechanism contribution 

has the form similar to that of MW-induced contribution 
to conductivity,® as it contains the factors sin(27ra;/u;c) 
and sin^(7ru;/wc). The first of these factors has extrema 
at oj/ojc = n ± 1/4, corresponding to the conditions of 
maximal displacement of electrons along the effective 
drag force or against this force under photon-assisted 
scattering, similar as in the case of a response to dc 
field.^ The second factor describes the enhancement of 
photon-assisted scattering probabilities in the resonance, 
oj/oJc = n, and their suppression in the anti-resonance, 
oj/oJc = n+ 1/2. The displacement-mechanism contribu¬ 
tion depend on MW polarization direction through the 
terms with the matrices of Eq. (38). 



9 


The fundamental difference between and the 

MW-induced contribution to conductivity is given by 
the factors Tj and Tj, which are not merely constants 
but functions of the magnetic field describing the mag¬ 
netophonon oscillations. The products of these magne¬ 
tophonon oscillating factors by the MW-induced oscillat¬ 
ing factors sin(27rw/u;c) and sin^(7ra;/u;c) physically cor¬ 
respond to the interference of these two kinds of oscilla¬ 
tions and can be viewed as a result of photon and phonon 
frequency mixing in the scattering probabilities. 

The phonon-drag part of /3 depends on electron tem¬ 
perature Te through the inelastic scattering time Tin ■ The 
quantities Fj and Fj are determined by the lattice tem¬ 
perature T. 

There is an important question of whether the ther¬ 
moelectric tensor /3 satisfies the symmetry with respect 
to time inversion (Onsager symmetry). In the absence of 
microwaves, this symmetry, of course, is satisfied. Under 
MW irradiation, when electrons are out of equilibrium, 
the Onsager symmetry can be broken. In application to 
the problem of electrons in the presence of electromag¬ 
netic waves, the time inversion implies, apart from the 
magnetic field reversal Wc —>■ — Wc, the transformations 
e — >■ e* and k —>■ —k, where k is the wave vector of the 
electromagnetic wave, e —>■ e* means that e± —>■ e^, 
which is equivalent to k± —>■ (see the beginning of 

Sec. II), while k —>■ — k means that the sign at tOp in 
Eq. (4) for s± is inverted, as follows from reversibility 
of the wave transmission problem^*^ employed for deriva¬ 
tion of Eq. (3). Therefore, the denominator in Eq. (4) 
transforms as a; ± Wc + —>■ a; =F Wc — which re¬ 

sults in s± —>■ s)p under the time inversion. The Onsager 
symmetry relation takes the form 

/3y(wc,s_,s+) =/3j*(-Wc,s%sl), (40) 

and similar relations can be written for the other trans¬ 
port coefficients including the conductivity. From s± —>■ 
one can see that both |s+p -I- |s_p and s_s!)_ are in¬ 
variants with respect to time inversion, thus the function 
b defined by Eq. (23) is also an invariant. The MW- 
induced part of P given by Eq. (37) does contain the 
terms violating the Onsager symmetry Eq. (40), these 
are the terms at the matrices go and gi. These terms 
change their signs under magnetic field reversal but are 
invariant under permutation of Cartesian indices. 


IV. THERMOPOWER COEFFICIENTS 

Having found P, one may calculate the thermopower 
tensor a, which is presented below as a sum of dark and 
MW-induced parts, a = . Because of the 

presence of terms which depend on MW polarization, this 
tensor is a general matrix. 

In the absence of MW irradiation, d has the same 
symmetries as the resistivity tensor, a^xx = ct^y and 

(0) (0) -O • • (0) 12 

ctxy = —o-yx- By using the expressions pxy = muicle. ns 


and pSx = mvtri^ + 2d'^)le'^ns together with Eqs. (31) 
and (32), where the Shubnikov-de Haas terms are ne¬ 
glected, one obtains, within the accuracy up to dP, the 
following results: 


= - 


3|e| ep 




-^(Fi + 2d2r,i) 

e 


(41) 


^ (0) _ ^c^tr 

~ ,.,2 I 


3|e| £f 


2d^ Ztx - 


2ttZ 


Wc T 


- ^2d"F,i 
e 


(42) 


In the classical case, the thermopower coefficients have 
usual forms found in literature.The Landau quanti¬ 
zation leads to additional terms proportional to d^. In 
the phonon-drag part of thermopower, these terms are 
determined by the function Fd oscillating with the mag¬ 
netic field. Because of these quantum corrections, the 
transverse phonon-drag thermopower is nonzero. In Fig. 
1 the longitudinal and transverse thermopower are plot¬ 
ted as functions of magnetic field for a rectangular GaAs 
quantum well of width 14 nm, with electron density 
Us = 5 X 10^^ cm“^ and mobility 2 x 10® cm^/V s. The 
quantum lifetime r = 7 ps is assumed, which corresponds 
to the ratio Tj^/t ~ 11. The phonon scattering time t\ 
is chosen as 0.2 gs for each mode, which approximately 
corresponds to 1 mm mean free path for phonons.®® The 
elastic coefficients for GaAs in units 10^^ dyn/cm^ are 
cii = 12.17, Ci 2 = 5.46, and C 44 = 6.16. The deforma¬ 
tion potential, piezoelectric coefficient, and density are 
V = 7.17 eV, hi 4 = 1.2 V/nm, and p = 5.317 g/cm®, re¬ 
spectively. The energy dependence of the transport time 
and quantum lifetime is assumed to be oc e®/^ and oc 
respectively, which corresponds to v{q) oc exp(—Zc?) un¬ 
der the condition of small-angle scattering, when IcPf ^ 
1. The oscillations of the thermopower coefficients are 
caused by magnetophonon resonances. At low tempera¬ 
ture, the oscillations are barely visible because the sys¬ 
tem falls into the Bloch-Gruneisen regime, but they are 
essential at higher temperatures. The last peak of a^xx is 
due to the scattering of electrons by high-energy (longitu¬ 
dinal) phonons, this peak disappears first with lowering 
temperature. The non-oscillating, proportional to 1/B, 
part of a^xy is determined by the diffusive contribution. 

Let us consider now the thermopower coefficients in the 
presence of MW excitation. While Eqs. (41) and (42) are 
valid for both classically strong and classically weak mag¬ 
netic fields, the MW-induced contributions are important 
only in the limit of classically strong magnetic fields. For 
this reason, only a part of the terms presented in Eqs. 
(35) and (37) are essential for calculation of thermopower 
in this limit. In particular, the longitudinal thermopower 
in classically strong magnetic fields is written simply as 
cXxx = Pxy Pyx- The influence of microwaves on Hall resis¬ 
tivity Pxy is weak^^, so Uxx is directly determined by Pyx- 
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B (T) B (T) 


FIG. 1: (Color online) Longitudinal (a) and transverse (b) 
thermopower in the absence of microwave excitation plotted 
for three temperatures. The calculations are done for a GaAs 
quantum well of width 14 nm with electron density ris = 
5 X 10^^ cm“^ and mobility 2 x 10® cm^/V s. The quantum 
lifetime of electrons is r = 7 ps. 


Neglecting the contributions of higher order in vtrl^c in 
Eqs. (35) and (37), one obtains 


(MW) _ 


2<PPu, 


2i:uj . 




■ Sin ■ 


X I r^i - 


OJc OJc 

'Tin^c \ 


STT'^epTeJ 


+ [(1 + b')Tc 2 + b'Tc 2 ] sin^- 

UJc 

+ [(1 + b')rs2 + 6'fs2]—sin^^l , 

UJc UJc j 


(43) 


while differs from this expression by changing 

the sign at b'. The first term in Eq. (43) is caused 
by modification of the isotropic distribution function of 
electrons by microwaves (inelastic mechanism) and in¬ 
cludes both the phonon-drag and the diffusive contribu¬ 
tions. Since the diffusive term increases with decreasing 
temperature, it may become comparable to the phonon- 
drag one. However, inevitable heating of electron gas by 
microwaves tends to hinder the contribution of the diffu¬ 
sive term. The remaining terms in Eq. (43) describe the 
phonon-drag thermopower caused by the displacement 
mechanism. They contain contributions proportional to 
6', which change the symmetry of the thermopower co¬ 
efficients. The dependence of these contributions on the 
polarization angle x can be illustrated for the case of 
linear polarization of the incident wave, when b is repre¬ 
sented in the form 


b = 


1 _ 2 jx ‘2‘iuJcUJp 

2® -b a;2 -b ujI 


(44) 


Since b' = Re(5) contains the terms both even and odd 


in magnetic field, Uxx, in general, is not symmetric in 
B (the reversal of magnetic field means alteration of 
the sign of uJc in all equations). The ’’inelastic” con¬ 
tribution in Eq. (43) should dominate at low enough 
temperatures, when vtrTin > 1- The ’’displacement” 
terms become more important with increasing temper¬ 
ature. It is worth to emphasize that the oscillations in 
these terms due to the factor sin^(7rw/a;c) are compara¬ 
ble by amplitude with the oscillations due to the factor 
sin(27ra;/a;c). This behavior is in contrast with that for 
MW-induced resistance. In the resistance, the contri¬ 
bution at sin(27ra;/a;c) dominates because it overcomes 
the oscillating part of sin^(7rw/wc) by the factor ^Txjjjtjjc 
which is numerically large in the region uj > ujc where 
MIRO are observed. As a consequence, the MW-induced 
resistance magnetooscillations due to the displacement 
mechanism are very similar to the magnetooscillations 
due to the inelastic mechanism,® so these two mechanisms 
are difficult to separate experimentally. In the phonon- 
drag thermopower, the contributions at sin(27ra;/a;c) and 
sin®(7ra;/a;c) are proportional to the functions rs 2 and 
rc 2 , respectively, and rc 2 is larger than rs 2 . More¬ 
over, rc 2 ^ rs 2 in the region of low magnetic fields, 
\ujc\ 47rpi?SA, see the Appendix. The ratio of the am¬ 
plitudes of sin(27ra;/a;c) and sin® (ttw/wc) oscillations in 
the ’’displacement” part of the thermopower is estimated 
as uj/2pps\, which is of the order of unity for typical 
electron densities and MW frequencies. The same is true 
for the ’’displacement” contribution to transverse ther¬ 
mopower described below by Eq. (47). 

A more careful analysis is required for evaluation of 
the transverse (Nernst-Ettingshausen) thermopower, be¬ 
cause the latter is determined by both diagonal and non¬ 
diagonal parts of /3 and is sensitive to MW-induced mod¬ 
ifications of the longitudinal resistivity. Indeed, axy = 
PxyPyy + PxxPxy The influence of microwaves on pxy is 
weak and not essential for determination of axy, while 
their influence on pxx is strong. Under the assumed con¬ 
dition that the electron-impurity scattering is more im¬ 
portant than electron-phonon scattering, the longitudinal 
resistivity correction due to MW irradiation is written as® 


JMW) ^ 
rxx 


2cPmpl 


e®ns 


trBn n 27ra; . 27ra; 

- Pu -sm-, 


UJc 


(45) 


and . Equation (45) implies that 

is governed by the inelastic mechanism. The displace¬ 
ment mechanism for electron-impurity scattering is less 
important at low temperatures, especially in the case 
of small-angle scattering processes relevant for high- 
mobility 2D systems.® In contrast, for electron-phonon 
scattering determining phonon-drag thermopower, the 
displacement mechanism is significant under the con¬ 
dition uj\Q < 2T when the main contribution to os¬ 
cillating functions rc 2 and rs 2 comes from large-angle 
scattering processes (backscattering). Among the ’’dis¬ 
placement” terms contributing into the transverse ther¬ 
mopower a^x!y^'^ there is a strong polarization-dependent 
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term coming from the diagonal part of the matrices 
go and gi in Eq. (37). The other contributions to 
contain a small factor vtrluJc- Out of them, only 
the ’’inelastic” ones can compete with the mentioned 
polarization-dependent contribution. Therefore, with the 
assumed accuracy up to <P , the result is written as a sum 
of two terms: 


sm(2x + ?7 b) + 




(46) 


where 


2^2 P 

Aaxy = ^|6| 


(rc2 - fc2)sin' 


7TUJ 


2 _ 

CO, 

2i:co 


+(rs 2 - fs 2 )— sin 

UJr 


(47) 


and 


OfP p ,y2 
in _ zn 
^xy I I 


eco. 


27TCO 27 T ( jJ 
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(- 


uj n~in 


2Te€F \27ra; a;, 


(jJc CJ, 

2iico 
sm ■ 


ri + 


— (1 — Ztr) COS 


isp 2 
2'kuj\ 


OJc 


(48) 


To obtain ay^^\ one should change the sign at the 
second term in Eq. (46). Since the effects under con¬ 
sideration are linear in MW intensity, the polarization- 
dependent term is a harmonic function of the doubled 
polarization angle; a similar angular dependence is ex¬ 
pected for electrical resistivity.^’^ This term is charac¬ 
terized by the amplitude Aaxy and the phase angle rye 
which are, respectively, a symmetric and an antisymmet¬ 
ric function of the magnetic field. For linear polariza¬ 
tion, when Eq. (44) is valid, the phase angle is defined as 
tan rye = 2wcWp/(a;^ — uj'^ + ojp). One may introduce the 
effective polarization angle xb = X + ’7b/ 2 describing the 
direction of the ac electric field in the 2D plane, which is 
different from the polarization of the incident wave. The 
polarization-dependent term, in general, is not antisym¬ 
metric under reversal of B, though for special orientation 
of the incident ac field along x or y axes the symmetry 
property is preserved. If the 

angle xb is equal to tt/ 2 or 0, which means that the elec¬ 
tric field in the 2D plane is polarized along y or x axes (i.e. 
along or perpendicular to the temperature gradient), the 
polarization-dependent term is equal to zero. The contri¬ 
bution of this term can be experimentally distinguished 
from the other contributions by its dependence on the 
polarization. 

The polarization-independent term given by Eq. (48) 
contains several contributions of different origin, though 
all of them are caused by the inelastic mechanism. The 
first part [the second line of Eq. (48)] comprises three dif¬ 
ferent contributions. The first one, at Ti, comes from the 
MW-induced correction to resistance if the thermoelec¬ 
tric current is due to the phonon-drag mechanism. The 


second contribution comes from the MW-induced correc¬ 
tion to resistance if the thermoelectric current is due to 
the diffusive mechanism. These two contributions can 
be distinguished from each other by their temperature 
dependence. At low temperatures (roughly estimated as 
Te < 0.5 K), the second contribution can exceed the hrst 
one, as it decreases with slower [see Eq. (All) for low- 
temperature behavior of Tij. However, the MW heating 
of electron gas renders this regime practically unrealiz¬ 
able. The third contribution, at Tgi, is caused by the 
MW-induced correction to the phonon-drag part of ther¬ 
moelectric tensor. In contrast to the first and second 
contributions, this one contains magnetophonon oscilla¬ 
tions. However, in the region of helds where these oscilla¬ 
tions exist, [wcl < 2pps\, the term Tsi/2 is much smaller 
than r 1 . The second part [the last line of Eq. (48)] con¬ 
tains the contributions due to MW-induced correction to 
diffusive part of thermoelectric tensor. This part does 
not exceed the contribution proportional to FTe/isp in 
the second line of Eq. (48) under the assumed comdi- 
tion [wcl <C 27r^Te. Therefore, the contribution propor¬ 
tional to Ti dominates over the others in Eq. (48) in 
the relevant region of parameters. This means that mag¬ 
netooscillations of a™ are determined only by the ratio 
oj/lOc and are similar to MIRO. The magnetoocillations of 
the polarization-dependent term are more complicated, 
because they also have the magnetophonon constituent 
due to the factors rc 2 — rc 2 and rs 2 — rs 2 , [see Eq. (47), 
Fig. 6 and its discussion below]. Therefore, the two terms 
in Eq. (46) can be distinguished from each other not only 
by polarization dependence and H-inversion symmetry 
but also by the behavior of magnetooscillations. 

It is important to emphasize that the components of 
the thermopower tensor given by Eqs. (43) and (46) do 
not violate the Onsager symmetry. This fact requires an 
explanation in view of the observation (see the end of Sec. 
HI) that some terms in /3 violate this symmetry. Indeed, 
a is formed as a result of matrix multiplication of p and 
(3 and its full form does contain the terms violating the 
Onsager symmetry. However, such terms are small in 
comparison to the terms included in Eqs. (43) and (46), 
so they are neglected. 

Coming to presentation of numerical results, let us con¬ 
sider first the diffusive contribution to thermopower co¬ 
efficients. This contribution is given by Eqs. (41), (42), 
(43), and (46), where all T^ and E^ are set to zero. The 
inelastic scattering time here and below is estimated ac¬ 
cording to® Tin = The diffusive thermopower is 

not sensitive to MW polarization. The longitudinal dif¬ 
fusive thermopower axx is modified by the microwaves 
in two ways: through the heating of 2D electrons and 
through the quantum correction in Eq. (43). The calcu¬ 
lations (see Fig. 2) demonstrate that the heating mech¬ 
anism is more essential. In particular, it leads to a peak 
at cyclotron absorption frequency and to oscillations at 
small B caused by the oscillations of absorbed MW 
power due to Landau quantization. The transverse dif¬ 
fusive thermopower axy, in contrast, is considerably af- 
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fected by the MW-induced quantum corrections from Eq. 
(48). Among these corrections there is a term 
whose oscillations directly reproduce the MIRO pattern 
shown in the inset of Fig. 2. The calculations demon¬ 
strate that the other terms, those in the last line of Eq. 
(48), are equally important, although their contribution 
becomes weaker with increasing temperature. 



B (T) B (T) 



FIG. 3: (Color online) Microwave-induced corrections to lon¬ 
gitudinal thermopower at T = 1.5 K and T = 4.2 K due to 
inelastic (a) and displacement (b) mechainsms, for linearly 
polarized MW excitation of frequency 130 GHz and electric 
field Auj = 2 V/cm. The parameters of the system are the 
same as in Fig. 1. Two plots for T = 4.2 K in (b) correspond 
to two angles of MW polarization. 


FIG. 2: (Color online) Longitudinal (left) and transverse 
(right) diffusive thermopower at T = 1.5 K and T = 4.2 
K under the linearly polarized MW excitation of frequency 
130 GHz and electric field = 2 V/cm. The parameters of 
the system are the same as in Fig. 1. The dashed lines show 
the dark thermopower (no MW excitation). The narrow solid 
line in the right-hand part shows the result of approximation 
for T = 1.5 K. The inset presents the 
calculated behavior of the longitudinal resistance. 

Consider now the influence of microwaves on the ther¬ 
mopower coefficients in the presence of both diffusive and 
phonon drag mechanisms. Theoretical and experimental 
studies of GaAs quantum wells show that for the temper¬ 
atures above 0.5 K the phonon-drag contribution domi¬ 
nates over the diffusive one. Consequently, the behavior 
of thermopower is governed mostly by the influence of 
MW excitation on the phonon-drag contribution. For 
the typical parameters of MW excitation, the oscillat¬ 
ing quantum corrections given by Eq. (43) are of the 
order of several /iV/K. The partial contributions due to 
inelastic mechanism (the first term in Eq. (43)) and dis¬ 
placement mechanism (the remaining terms) are shown 
in Fig. 3. The role of the displacement mechanism in¬ 
creases with increasing temperature. At low tempera¬ 
tures (Bloch-Gruneisen regime), the period of the oscil¬ 
lations is determined by the ratio w/wc- With increas¬ 
ing temperature, the magnetophonon resonances become 
important and the picture of oscillations becomes more 
rich. The sensitivity of the displacement mechanism to 
MW polarization is illustrated by plotting its contribu¬ 
tion for two angles of electric field of the incident wave, 
X = 0 and x = 7r/4. 

However, the realtive change of the longitudinal com¬ 


ponent axx under MW irradiation is not strong. The 
terms due to phonon drag in Eq. (43) are proportional 
to the functions Tsi, rc 2 , and rs 2 , which are small in 
comparison to Ti in the important region of parame¬ 
ters |wc| < “^pfSxq and Iwd ^ir'^Te, where magne¬ 
tophonon oscillations take place but Shubnikov-de Haas 
oscillations are suppressed (see a more detailed compar¬ 
ison in the Appendix). The ratio of the relative change 
of axx due to MW irradiation to the relative change of 
the resistivity pxx is estimated by a small factor Tsi/Ti. 
This means that even in the case when MW-induced re¬ 
sistance oscillations are strong, the MW-induced oscilla¬ 
tions of the longitudinal thermopower still may be weak. 
The magnetic-field dependence of axx at low tempera¬ 
ture is presented in Fig. 4 (a). For T = 1.5 K one can 
see changes in the oscillation picture, in particular, inver¬ 
sion of the minimum around 0.18 T and a considerable 
enhancement of the last peak. The vertical shift of axx 
as a whole with respect to a^xx is caused by the diffusive 
mechanism contribution, due to heating of electrons by 
microwaves, see Fig. 2. With increasing temperature, 
the relative effect of microwaves on axx becomes weaker 
because ai^J increases faster than aiif. 

The transverse thermopower axy, in contrast, is 
strongly changed by microwaves, because the dark ther¬ 
mopower ai^J is small itself. At low temperature [see 
Fig. 4 (b)] the modification is almost entirely governed 
by the oscillations of resistivity, which means that the 
approximation ~ works well. This 

approximation is no longer valid when temperature in¬ 
creases and the polarization-dependent contribution, the 
first term in the expression Eq. (46), becomes significant. 
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FIG. 4: (Color online) Longitudinal (a) and transverse (b) 
thermopower at T = 1.5 K under the linearly polarized MW 
excitation of frequency 130 GHz and electric field E^, = 2 
V/cm. The parameters of the system are the same as in 
Fig. 1. The dashed lines show the dark thermopower. The 
narrow solid line shows the result of approximation ~ 

for transverse thermopower. 


This is demonstrated in Fig. 5, where axy is plotted for 
two directions of ac electric field: along x axis (x = 0) 
and at the angle of 7r/4 to this axis. With increasing B, 
when the ratio Vtrjojc. becomes smaller, a^y deviates from 

the simple dependence oc and becomes strongly 

sensitive to polarization. 



B(T) 


FIG. 5: (Color online) Transverse thermopower at T — 4.2 
K under the MW excitation of frequency 130 GHz and elec¬ 
tric field Eu, = 2 V/cm, for two different directions of linear 
polarization of incident wave. The parameters of the system 
are the same as in Fig. 1. The dashed line shows the dark 
thermopower. The narrow solid line shows the result of ap- 
proximation aiy ~ Pxx Pxy ■ 


The polarization dependence of a^y for different mag¬ 
netic fields is characterized by the amplitude Aaxy given 
by Eq. (47). This function is plotted in Fig. 6 for 
different temperatures. The complicated oscillating be¬ 
havior of Aaxy is caused by the interference of mag¬ 
netophonon oscillations with microwave-induced oscilla¬ 
tions. At small T, when the system is in the Bloch- 
Gruneisen regime, Aaxy is small. With increasing T, 
Aaxy increases and saturates around 10-15 K. The inset 
in Fig. 6 shows how the rotation of the MW polarization 
angle changes the total transverse thermopower. 

The relative contribution of polarization-dependent 
part can be further enhanced at higher MW intensity 
and at higher mobility, because the second term in Eq. 
(46) is proportional to the factor VirEn which goes down 
when inelastic scattering time oc T~^ decreases be¬ 
cause of microwave heating of electron gas and when the 
transport scattering rate vtr (inversely proportional to 
the mobility) decreases. 



B(T) 


FIG. 6: (Golor online) Magnetic-field dependence of 

polarization-sensitive part of transverse thermopower at dif¬ 
ferent temperatures, for the MW excitation of frequency 130 
GHz and electric field E^x = 2 V/cm. The parameters of the 
system are the same as in Fig. 1. The inset shows dependence 
of thermopower on the polarization angle at B = 0.5 T. 

In the case of circular polarization or non-polarized ra¬ 
diation (chaotic polarization) the polarization-dependent 
term vanishes and a^^^ is determined by the second 
term in Eq. (46). Since the most important part of this 

term is given by the oscillations of transverse 

thermopower under these conditions follow the MW- 
induced resistance oscillations. 

The longitudinal and transverse thermopower compo¬ 
nents axx and axy are directly measured in the Hall bars. 
The longitudinal thermopower can also be measured in 
the Corbino disc geometry.^® In this case, polarization- 
dependent terms do not appear and the voltage be¬ 
tween inner and outer contacts is determined by the ther- 
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mopower ad = Pd/crd, where Pd and ad are the diagonal 
parts of the tensors $ and a in the absence of MW po¬ 
larization. Since ad is modified by microwaves stronger 
than Pd, the behavior of thermopower in MW-irradiated 
Corbino discs is determined mostly by MW-induced os¬ 
cillations of ad- 

The theory developed in this paper does not take into 
account temperature dependence of the density of states. 
Such a dependence appears mostly due to contribution of 
electron-electron scattering into the inverse quantum life¬ 
time 1/t (see Ref. 11 and references therein). This effect 
leads to an exponential suppression of all quantum con¬ 
tributions in the transport coefficients, including those 
considered above, with increasing Tg. Formally, this oc¬ 
curs because the Dingle factor d acquires a multiplier 
exp(-7r/Tee(Te)|a;c|), where l/Tee{Te) ^ T^/ep- This ef¬ 
fect tends to decrease the quantum part of dark ther¬ 
mopower and MW-induced corrections to thermopower 
with increasing temperature. Since the main (phonon- 
drag) contribution to thermopower, in contrast, increases 
with increasing temperature at T < pfS\, it is impor¬ 
tant to investigate possible competition of these oppo¬ 
site trends in the quantum (proportional to dP) terms 
in thermopower. Assuming that Tg — T, the expo¬ 
nential dependence of these terms on temperature in 
the Bloch-Gruneisen regime {T <C PfS\) is written as 
where <I)t — 2ttT‘^/efIuJc] + 2pFS\/T is a non¬ 
monotonic function of temperature. This function de¬ 
creases at T < Tq and increases at T > Tq, where 
To ~ pFSA(|wc|/47rTOs|)^/^. Since the estimate for GaAs 
gives To > pfS\ even for magnetic fields as small as 0.05 
T, one may conclude that the temperature dependence 
of the density of states does not alter the thermal in¬ 
crease of the quantum contributions to thermopower at 
T < pfS\. However, at T > pfS\ all these contribu¬ 
tions, both in the dark thermopower and MW-induced 
corrections, decrease with temperature instead of going 
to saturation. 


V. DISCUSSION AND CONCLUSIONS 

The influence of MW irradiation on the energy distri¬ 
bution of electrons and on electron scattering by phonons 
and impurities has a profound effect on transport prop¬ 
erties of 2D electron systems in perpendicular magnetic 
field. While the effect of microwaves on the electrical 
resistance is widely studied, the related behavior of the 
other kinetic coefficients has not received proper atten¬ 
tion. This paper reports a theoretical study of possible 
MW-induced quantum effects in thermopower. Such ef¬ 
fects can exist in the samples with high electron mobility 
in the moderately strong magnetic fields, that is, under 
the same conditions when the MW-induced quantum os¬ 
cillations of the electrical resistance are observed. 

In contrast to electrical resistance, which at low tem¬ 
peratures is determined by electron-impurity scatter¬ 
ing, the thermopower is determined mostly by electron- 


phonon scattering, through the phonon drag mechanism. 
The theory of phonon-drag thermoelectric response in 
quantizing magnetic fields remains an issue of interest 
even under quasi-equilibrium conditions, in the absence 
of MW irradiation. A further development of such theory 
is presented in this paper. In particular, an anisotropy 
of the acoustic phonon spectrum has been taken into ac¬ 
count and analytical expressions valid in the regime of 
overlapping Landau levels with the accuracy up to the 
square of the Dingle factor have been derived, see Eqs. 
(32), (33), (41) and (42). The theory gives a clear picture 
of the origin of magnetophonon oscillations observed^® 
in the longitudinal thermopower of high-mobility GaAs 
quantum wells and predicts similar oscillations in the 
transverse thermopower (Fig. 1). For typical param¬ 
eters of GaAs wells, the oscillations are clearly visible 
for temperatures above 2 K, while at lower temperatures 
they become exponentially suppressed because the Bloch- 
Gruneisen regime is reached. In the experiment,how¬ 
ever, the oscillations were resolved between 0.5 K and 
1 K. This discrepancy can be explained by taking into 
account that the phonon distribution function in the ex¬ 
periments on thermopower is not reduced to the form of 
Eq. (17) commonly applied by theorists. Even at low 
temperatures of the sample, there can exist high-energy 
phonons able to cause backscattering of electrons. In¬ 
deed, since the phonon mean free path at low tempera¬ 
tures is very large (of 1 mm scale), it is quite possible 
that such high-energy phonons may arrive to the 2D sys¬ 
tem directly from the heater, via ballistic propagation. 
Another possible reason, which is especially relevant at 
low temperatures, is that the modification of phonon dis¬ 
tribution function is strong and cannot be represented 
in the form of a small correction linear in temperature 
gradient. In any case, a quantitative agreement with ex¬ 
periment can be reached only if the phonon distribution 
is known. The theory presented in this paper can be 
generalized to the case of arbitrary phonon distribution 
by substituting the antisymmetric part of actual phonon 
distribution function instead of the second term in Eq. 
(17). 

The influence of MW irradiation on the longitudinal 
axx and transverse axy components of the thermopower 
has been studied above by using the approved meth¬ 
ods applied earlier to calculation of the resistivity. It 
is found that the MW irradiation has a considerable ef¬ 
fect on both these components. In contrast, for elec¬ 
trical resistance the microwaves strongly modify only 
the longitudinal component pxx- Both the diffusive and 
phonon-drag contributions to thermopower are shown to 
be affected by MW irradiation. The MW-induced quan¬ 
tum corrections to diffusive thermopower increase with 
decreasing electron temperature, in contrast to classi¬ 
cal diffusive thermopower, which is proportional to this 
temperature. However, since the phonon-drag contribu¬ 
tion dominates, the MW-induced quantum corrections to 
phonon-drag thermopower appear to be more important. 
These effects are of the order of several pV /K for typical 
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parameters of the 2D system and MW excitation, and 
can be detected experimentally. The oscillating behav¬ 
ior of MW-induced corrections as functions of the mag¬ 
netic field reflects the properties of electron scattering 
by phonons under conditions when the electron distri¬ 
bution function acquires MW-induced oscillating compo¬ 
nent (inelastic mechanism) and when MW-assisted scat¬ 
tering takes place (displacement mechanism). Both these 
mechanisms are important, and both provide a mixing 
of resonant phonon frequencies with MW frequency w, 
thereby leading to interference oscillations of the ther¬ 
mopower. 

In terms of relative values, the MW-induced changes 
in the longitudinal thermopower are much smaller than 
the corresponding effect in the resistivity. In contrast, 
the relative MW-induced changes in the transverse ther¬ 
mopower are large, because in the classically strong mag¬ 
netic fields the transverse thermopower itself is much 
smaller than the longitudinal one. At lower temperatures 
and weaker magnetic fields, the oscillations of transverse 
thermopower a.xy follow the picture of MW-induced re¬ 
sistance oscillations (MIRO) [Fig. 4 (b)]. As the temper¬ 
ature and magnetic field increase, the oscillations of o-xy 
no longer follow the MIRO picture and become strongly 
sensitive to polarization of the incident wave. The po¬ 
larization dependence of Uxy is much stronger than the 
corresponding dependence of the electrical resistivity un¬ 
der MW irradiation. These finding may stimulate exper¬ 
imental studies of the transverse thermopower of MW- 
irradiated 2D electron gas. 

The appearance of a large polarization-dependent term 
in the MW-induced transverse thermopower is one of the 
main results of the present study. The nature of this 
effect can be easily understood by considering the colli¬ 
sionless approximation (no electron-impurity scattering, 
vtr = 0), when the transverse thermopower does not ap¬ 
pear without MW irradiation. The drag of electrons by 
the phonons drifting along the temperature gradient VT 
can be described^^ in terms of a dragging force due to 
effective electric field Ep/j oc VT. The electrons in the 
magnetic field are drifting perpendicular to Ep^,. To com¬ 
pensate this drift, a real electric field E = —Ep?i develops. 
Thus, the longitudinal thermopower is equal to |E|/|VT| 
while the transverse thermopower is zero. When a polar¬ 
ized ac field is applied to the system, the effective elec¬ 
tric field Ep/i, in general, is not directed along VT and 
becomes sensitive to polarization. This occurs because 
Ep/i is formed as a result of electron-phonon interaction 
assisted by emission and absorption of radiation quanta, 
and this interaction is stronger when the in-plane compo¬ 
nents of phonon momenta are parallel to the polarization- 
dependent vector see Eqs. (12) and (13). Conse¬ 
quently, the real electric field E = — Ep?i is not paral¬ 
lel to VT, which means that there exists a transverse 
component of thermopower. This component is given by 
the first term in Eq. (46). Beyond the collisionless ap¬ 
proximation, the other, polarization-independent terms 
in (Xxy are also important. A larger relative contribution 


of polarization-dependent term is expected in 2D electron 
systems with higher mobility (smaller vtr)- 

An important issue left beyond the above considera¬ 
tion is the behavior of thermopower at zero longitudinal 
resistance. In high-mobility 2D systems, intensive MW 
irradiation leads to a remarkable phenomenon of zero 
resistance states,which means that the longitudinal 
resistance vanishes in certain intervals of magnetic 
fields corresponding to MIRO minima at lower MW 
intensity. This effect is often explained (see Ref. I 
and references therein) as a result of the instability of 
homogeneous current flow under condition of negative 
local resistance, which leads to spontaneous formation 
of domains with different directions of the currents and 
Hall fields. Since the longitudinal resistivity formally 
enters the expression for thermopower and, as shown 
above, considerably affects the transverse thermpower 
in the presence of MW irradiation, the magnetic-field 
dependence should demonstrate the regions of nearly 
constant ocxy in the intervals of pxx = 0, while o-xx is 
not expected to be sensitive to zero resistance states. Of 
course, this conclusion looks somewhat naive, because 
the presence of domains may affect the behavior of mea¬ 
sured thermopower. It is not clear, however, which kind 
of domain picture is realized under zero resistance state 
conditions in thermoelectric experiments, when there is 
no electric currents through the contacts. Future studies 
should shed light on this particularly interesting problem. 
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Appendix A: Asymptotic behavior of the functions 

Ti and Ti 


In the approximation of isotropic phonon spectrum, 
the integral over the polar angle (pq in the operator Vn 
can be carried out analytically, and Eq. (33) is reduced 
to the form 



dO 

Pm Jo 



dqz 

TT 


X ^ rxGxF 



cos 

OJc 

27TSxQ 



sin 


(Al) 


where Q = + q = 2pFSm{0/2), Gi = 

+ {ehiiY9(Gql/2Q^, and Gt = (e^i4)^(8g^gf -I- 
q^)/2Q^. Eor fi one should replace Gi and Gt by G; = 
-{ehufdq'^ql/4:Q^ and Gt = {ehuf {Sq'^ql - q^)/AQ^, 
respectively. The functions Gt and Gt describe interac¬ 
tion of electrons with transverse phonon modes due to 
piezoelectric potential mechanism, while G; and G; de¬ 
scribe interaction with longitudinal phonon modes due 
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to both deformation potential and piezoelectric poten¬ 
tial mechanisms. Analytical expressions for the functions 
hen; Ten, and Fgn calculated from Eq. (Al) are 
given below in some limiting cases. 

In the limit Wc <C ^'ksxPf, when cos{2'ks\Q/lOc) and 
s\n{2TTsxQ/ lOc) are rapidly oscillating functions of 9 and 
9zIpf, the main contribution to the integrals in Eq. (Al) 
comes from the region of small when Iq^ ~ 1, and 
from two regions of 6 around 0 = 0 (corresponding to 
forward scattering of electrons) and 9 = n (backscatter- 
ing), because these are the regions of most slow variation 
of (5 as a function of 9 and g^. Under the requirement 
\ojc\ ^ 2 tt‘^T, which is already stated as the condition 
when the Shubnikov-de Haas oscillations are suppressed, 
one obtains 
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Eor comparison, it is useful to present also the expres¬ 
sion for Fi: 


Fi 


177C(3) , ( 
29 5^* 



where (’(fc) is the Riemann zeta-function. This expression 
is valid in the limit of T <C sxPf and can be used for 
order-of-value estimates at T ~ sxPf- 

From the definition (AlO), the applicability region for 
Eqs. (A2) - (A9) can be written as ca ^ 1. The 
magnetooscillations of the functions described by Eqs. 
(A 2 ) - (A9) occur because of the terms with cosca and 
sincA. The amplitudes of these oscillating terms are al¬ 
ways much smaller than Fi of Eq. (All) in the case 
ca ^ 1. If T ~ saPF) this smallness is given by the 
factors for Fd, Fc 2 , Fci, and rc 2 and e~^'^ for F^i, 
Fs 2 , Fsi, and rs 2 . With lowering T, the oscillations 
are exponentially suppressed because of F{sxPf/T) ~ 
{2sxPf/T)"^ eyip{—2sxPF/T) at T <C sxPf- In the case 
of strong exponential suppression, the absolute values of 
the functions given by Eqs. (A2) - (A9) are determined 
by their non-oscillating parts which are proportional to 
powers of uje- The non-oscillating parts of n = 1 func¬ 
tions (Fci, Fgi, Fci, and F^i) are much smaller than Fi 
due to parameters (wc/27r^T)^ for piezoelectric-potential 
contribution and (wc/27r^T)^ for deformation-potential 
contribution. The non-oscillating parts of n = 2 func¬ 
tions (Fc 2 , Fs 2 , Fc 2 , and Fs 2 ) contain extra small factors 
e)C^, because these functions are much smaller than n = I 
functions at small-angle scattering, 0 ^ 1 . 

In stronger magnetic fields, when oje is comparable 
to 47rsA_PF, analytical expressions can be obtained at 
T > SxPf and under a wide-well approximation, the lat¬ 
ter means that the quantum well width a is much larger 
than F jpF so that the convergence of the integral over g^ 
takes place at gz <C pf and is governed by the function 
Iq^. Introducing go = dqzlq^ (for infinitely deep 

rectangular well go = 3/2a), one obtains 
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Fi is given by 
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Xk{x) 


d’^Jojx) 

dx^ 


(A16) 


is the k-ih order derivative of the Bessel function Jo{x). 
Such derivatives can be expressed through the other 
Bessel functions Ji{x). In the special case of Fsi, there is 
a term with the function which should be treated 

as the antiderivative of Jo(et)- This term is expressed 
through the Bessel functions and Struve functions H^: 
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In the regime of validity of Eqs. (AI2)-(AI5) the function 
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For large arguments the functions (A12)-(A15) are 
reduced to combinations of oscillating factors sine^ and 
cosca, similar to the case described by Eqs. (A2)-(A9), 
and are small in comparison to Fi. If — I, these func¬ 
tions become comparable to Fi. Actually, the wide well 
limit a 3> 'xjpp is hardly attainable for single-subband 
occupation in the quantum well. The expressions (AI2)- 
(A15) are nevertheless useful for estimates of the maxi¬ 
mal possible values of the quantities F^ and F^. 
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